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We propose a systematic, model-independent treatment of electric dipole transitions of heavy 
— 1 , quarkonium. Within an effective field theory, concretely potential non-relativistic QCD, the 

relativistic corrections of relative order v 2 to the decay rate are derived. An existing formalism 
(-H ' developed for Ml decays will be extended for our purpose. We scrutinize and complement 

C^H 1 former results from potential model calculations. 

& 

1 Introduction 

>. 

£\j | Radiative transitions play an important role for our understanding of QCD, in particular 

£-». ' of heavy quarkonium. They provide information about the wave functions describing the 

physical system and probe both the perturbative and non-perturbative regime. Especially 

El transitions give significant contributions to the total decay rate and are observed in the 

experimental facilities. Recently decays of charmonium were measured at BES and CLEO, 

t— I | including the observation of the process h c — >■ // c 7 in 2010 [1]. Concerning bottomonium, 

J> ■ CLEO, BaBar and Belle produced many data, e.g. for the determination of the branching 

fractions Xb states [2,3]. A review about recent developments can be found in [4]. 

On the theory side, electric dipole transitions were treated in several potential models, a 
summary can be found in [5]. We will refer to [6] for comparison with the general results 
for ^-decays. A model-independent treatment to check and improve the calculations has 
been missing so far. However, in the last decade there has been significant progress using 
effective field theories (EFTs) to describe heavy quarkonium (see [7] and references therein). 
Since heavy quarkonium is assumed to be a non-relativistic system we may take advantage 
of the hierarchy of scales m 3> mv 3> mv 2 , where v <C 1 is the heavy quark velocity, m 
is the heavy quark mass ("hard scale"), p ~ mv is the relative momentum of the bound 
state ("soft scale") E ~ mv 2 is the binding energy E ~ mv 2 ("ultrasoft scale"). The ultimate 
EFT living at the ultrasoft scale is potential non-relativistic QCD (pNRQCD). In 2005, for 
the first time radiative decays, concretely Ml transitions, were calculated in this theory [8]. 
Using the framework of that paper as a guideline we close the gap and compute the decay 
rates of the El processes n 3 Pj — > n' Si 7 and n x V\ — > n' So 7. The following is based on [9]. 
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2 The Lagrangians in NRQCD and pNRQCD 

By integrating out the hard scale m 3> Aqcd from the fundamental theory, QCD, in per- 
turbation theory (« s (m) <C 1) one obtains non-relativistic QCD (NRQCD) [10,11]. For the 
calculation of El transitions at NLO only the two-fermion Lagrangian £2-/ matters and 
the relevant part reads 



/ r> 2 F) 4 \ / r em r em \ 

(D &-f = ^ ('Do + ^ + ^3) * + ee Q f [^cr ■ B«" + .^^[D^B*]] </> + ex. . 

with zD = ia - gT s Ag - eeQA^ 1 , z'D = z'V + gT^A" + ee Q A em and t/> denoting a Pauli 
spinor for the heavy quark. The matching coefficients are found to be c" n = 1 + Cpj^ + 
£>(a 2 ) / cf"=24 m -l. 

For processes at the ultrasoft scale, NRQCD is not yet the appropriate theory, since there 
are still several scales entangled (p,E, Aqcd) and thus no homogeneous power counting 
can be established. Integrating out the soft scale mv we obtain a theory for ultrasoft modes, 
i.e. pNRQCD [12, 13]. The crucial step to disentangle the energy and momentum scale is 
the multipole expansion in the relative distance r. To be definite we will work in the weak- 
coupling regime, where p>E> Aqcd- The power counting reads 

(2) r ~ 1/mv, V r = d/dr ~ mv, V = 9/9R ~ mv 2 , E, B ~ (mv 2 ) 2 , E em , B em ~ k 2 . 

k 7 is the energy of the emitted photon, which scales like fc 7 ~ mv 2 for transitions between 
states with different principal quantum numbers. 

The pNRQCD-Lagrangian contributing at NLO in the decay rate, i.e. at order k^v° /m 2 , 
reads 



+ gV A (0 + r-ES + S + r-EO)} 
(3) + £ 7 pNRQCD + flight / 

where the covariant derivatives are given by iD$0 = idoO — g[T a A a Q , O] and /DO = /VO + 
g[T fl A",0] and the trace goes over the color and spin indices. The singlet potential V$ 
has been calculated perturbatively and non-perturbatively to order 1/m 2 ( [14-16], for 
more original references see [7]), we display the structure of the relevant potentials for 
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pNRQCD £ , * + B - , <|_ 



Figure 1 : Example for a reducible diagram, if the electromagnetic operator commutes 
with the gluonic ones. It does not contribute to the matching coefficient of a single opera- 
tor. 



computations at NLO in the decay rate, 

(4) vi W =vm W + ^W + ^W + ^W. 

m m z m z 

(5) vg>(r) =vV {r ) + \{V$\r),f} + V l^l}, 

(6) Vjg(r) =V$(r)L • S + y s ( 2 2) (r)S 2 + vg(r) [3(f ■ «n)(f • cr 2 ) - a x ■ cr 2 ] . 
The relevant part of £ 7P nrqcd f° r El transitions is 

£?JnrQCd = ee Q I d 3 r7rf.V rE S f r ■ E em S + V£ E + r • E em O + -Ly('" v ) 2 '" £ S + r • [(rV) 2 E raj ]S 

+ t-LyV-(rxB) s tr v xB «»i s 
4m i / j 

+ fJ_yV r .(rx(rV)B) s t| Vrvrx[(rV)Bem] | S 



+ J-y('' v )^ B [S + / £r].[(rV)B eM ]S 



4m 

(7) _;_Ly^-(Exvv)[ S + /£r ] . ^ e m x v ^ s | _ 

In fact more terms are allowed according to the symmetries of pNRQCD. However, we can 
show that their matching coefficients vanish. The matching is done by equating Green's 
functions in NRQCD and pNRQCD at the energy scale mv order by order in the inverse 
mass and r. The crucial argument for several operators is that diagrams in NRQCD which 
can be cast into a reducible structure also give reducible diagrams in pNRQCD. Therefore 
they have to be subtracted to obtain irreducible operators in pNRQCD and do not play a 
role in the matching procedure. An example is the diagram in Fig. 1, where the gluonic 
contribution can be factorized out yielding just a potential. Using this argument we can fix 
all of the Wilson coefficients in (7), so that the exact QCD results reproduce the ones from 
tree level calculations, namely 
(8) 

yr-E _ yr-E _ y{rV) 2 r-E _ yV-(rxB) _ y(rV)V r -(rxB) _ -y y(rV)a-B _ Q em ya-(ExV r ) _ c em 
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Figure 2: Color octet contributions to El transitions, the double line represents the inter- 
mediate octet state. 



3 Results 

With the help of the formalism developed in [8] we can describe the states in a quantum 
mechanical way using wave functions and compute the decay rate at NLO from the La- 
grangian (7). We obtain 
(9) 



r„ P ,^ Sl7 = I «„44«("i - »'°) ( » + * - 1| - £ + 2 ,„ 



fc 7 (c e F m 



K/ + i) 



where 



(10) 



I N = I drr N R n , (r)R nl (r) 

o 



R contains all of the wave-function corrections due to the higher-order potentials men- 
tioned in (4)-(6), the relativistic correction of the kinetic energy — p 4 /4m 3 , and higher-order 
Fock state contributions which are given by the diagrams in Fig. 2. In contrast to Ml tran- 
sitions the latter ones do not vanish for El decays. 

The expression (9) is also valid in the strongly coupled regime (without color-octet contri- 
butions in R), where p ~ ^qcd> since we made use of non-perturbative matching argu- 
ments and additional operators do not appear in this regime. 

Compared to the results with the potential model calculation in [6] we find an equivalence 
between (9) and the corresponding formula there at the given order. However, our definite 
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power counting allowed us to include all relativistic corrections systematically, in partic- 
ular the color-octet contributions in the weak-coupling regime and the one coming from 
Vy . Both were missing in former approaches. Furthermore we can show that the anoma- 
lous magnetic moment c e ™ — 1 ~ 0(a s (m)) is actually suppressed and does lead to large 
non-perturbative contributions. 

Without much effort one can extend the discussion to other processes like n l P\ — > n' So 7 
and n Si — > n Pjj. Finally based on these results a phenomenological analysis for char- 
monium and bottomonium decays can be performed. 
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